********* Table A20 *********
clear all
set more off 
capture log close
 
* main path to data
global path "~/Dropbox/Paraguay/1_data/7_replication_dataverse"
capture cd "$path"	

*tables = analysis tables
global tables "$path/tables"

* load rural collective action data
use "$path/data/landinv.dta", clear

xtset district_code year

// Generate indicator variables for years 
forval year = 1954/2007 {
    generate y`year' = (year == `year')
}

// Generate the alt_1989 variable and set it to 0
generate alt_1988 = 0

set matsize 10000
* count number of observations in the sample
count if inrange(year, 1983, 1995)
local obs_shown = r(N)


set more off
eststo clear

reg l_linv  c.betweenness_z  y1955#c.betweenness_z ///
			y1956#c.betweenness_z y1957#c.betweenness_z y1958#c.betweenness_z y1959#c.betweenness_z y1960#c.betweenness_z ///
			y1961#c.betweenness_z y1962#c.betweenness_z y1963#c.betweenness_z y1964#c.betweenness_z y1965#c.betweenness_z ///
			y1966#c.betweenness_z y1967#c.betweenness_z y1968#c.betweenness_z y1969#c.betweenness_z y1970#c.betweenness_z ///
			y1971#c.betweenness_z y1972#c.betweenness_z y1973#c.betweenness_z y1974#c.betweenness_z y1975#c.betweenness_z ///
			y1976#c.betweenness_z y1977#c.betweenness_z y1978#c.betweenness_z y1979#c.betweenness_z y1980#c.betweenness_z ///
			y1981#c.betweenness_z y1982#c.betweenness_z y1983#c.betweenness_z y1984#c.betweenness_z y1985#c.betweenness_z ///
			y1986#c.betweenness_z y1987#c.betweenness_z alt_1988 y1989#c.betweenness_z y1990#c.betweenness_z ///
			y1991#c.betweenness_z y1992#c.betweenness_z y1993#c.betweenness_z y1994#c.betweenness_z y1995#c.betweenness_z ///
			y1996#c.betweenness_z y1997#c.betweenness_z y1998#c.betweenness_z y1999#c.betweenness_z y2000#c.betweenness_z ///
			y2001#c.betweenness_z y2002#c.betweenness_z y2003#c.betweenness_z y2004#c.betweenness_z y2005#c.betweenness_z ///
			y2006#c.betweenness_z y2007#c.betweenness_z ///
			i.year i.district_code  i.year##i.depa ,  cluster(district_code)

eststo
estadd scalar obs_shown = `obs_shown'

reg linv_b  c.betweenness_z  y1955#c.betweenness_z ///
			y1956#c.betweenness_z y1957#c.betweenness_z y1958#c.betweenness_z y1959#c.betweenness_z y1960#c.betweenness_z ///
			y1961#c.betweenness_z y1962#c.betweenness_z y1963#c.betweenness_z y1964#c.betweenness_z y1965#c.betweenness_z ///
			y1966#c.betweenness_z y1967#c.betweenness_z y1968#c.betweenness_z y1969#c.betweenness_z y1970#c.betweenness_z ///
			y1971#c.betweenness_z y1972#c.betweenness_z y1973#c.betweenness_z y1974#c.betweenness_z y1975#c.betweenness_z ///
			y1976#c.betweenness_z y1977#c.betweenness_z y1978#c.betweenness_z y1979#c.betweenness_z y1980#c.betweenness_z ///
			y1981#c.betweenness_z y1982#c.betweenness_z y1983#c.betweenness_z y1984#c.betweenness_z y1985#c.betweenness_z ///
			y1986#c.betweenness_z y1987#c.betweenness_z alt_1988 y1989#c.betweenness_z y1990#c.betweenness_z ///
			y1991#c.betweenness_z y1992#c.betweenness_z y1993#c.betweenness_z y1994#c.betweenness_z y1995#c.betweenness_z ///
			y1996#c.betweenness_z y1997#c.betweenness_z y1998#c.betweenness_z y1999#c.betweenness_z y2000#c.betweenness_z ///
			y2001#c.betweenness_z y2002#c.betweenness_z y2003#c.betweenness_z y2004#c.betweenness_z y2005#c.betweenness_z ///
			y2006#c.betweenness_z y2007#c.betweenness_z ///
			i.year i.district_code  i.year##i.depa ,  cluster(district_code)
			
eststo
estadd scalar obs_shown = `obs_shown'

reg l_peasfam c.betweenness_z  y1955#c.betweenness_z ///
			y1956#c.betweenness_z y1957#c.betweenness_z y1958#c.betweenness_z y1959#c.betweenness_z y1960#c.betweenness_z ///
			y1961#c.betweenness_z y1962#c.betweenness_z y1963#c.betweenness_z y1964#c.betweenness_z y1965#c.betweenness_z ///
			y1966#c.betweenness_z y1967#c.betweenness_z y1968#c.betweenness_z y1969#c.betweenness_z y1970#c.betweenness_z ///
			y1971#c.betweenness_z y1972#c.betweenness_z y1973#c.betweenness_z y1974#c.betweenness_z y1975#c.betweenness_z ///
			y1976#c.betweenness_z y1977#c.betweenness_z y1978#c.betweenness_z y1979#c.betweenness_z y1980#c.betweenness_z ///
			y1981#c.betweenness_z y1982#c.betweenness_z y1983#c.betweenness_z y1984#c.betweenness_z y1985#c.betweenness_z ///
			y1986#c.betweenness_z y1987#c.betweenness_z alt_1988 y1989#c.betweenness_z y1990#c.betweenness_z ///
			y1991#c.betweenness_z y1992#c.betweenness_z y1993#c.betweenness_z y1994#c.betweenness_z y1995#c.betweenness_z ///
			y1996#c.betweenness_z y1997#c.betweenness_z y1998#c.betweenness_z y1999#c.betweenness_z y2000#c.betweenness_z ///
			y2001#c.betweenness_z y2002#c.betweenness_z y2003#c.betweenness_z y2004#c.betweenness_z y2005#c.betweenness_z ///
			y2006#c.betweenness_z y2007#c.betweenness_z ///
			i.year i.district_code  i.year##i.depa ,  cluster(district_code)
			
eststo
estadd scalar obs_shown = `obs_shown'

reg l_hectares   c.betweenness_z  y1955#c.betweenness_z ///
			y1956#c.betweenness_z y1957#c.betweenness_z y1958#c.betweenness_z y1959#c.betweenness_z y1960#c.betweenness_z ///
			y1961#c.betweenness_z y1962#c.betweenness_z y1963#c.betweenness_z y1964#c.betweenness_z y1965#c.betweenness_z ///
			y1966#c.betweenness_z y1967#c.betweenness_z y1968#c.betweenness_z y1969#c.betweenness_z y1970#c.betweenness_z ///
			y1971#c.betweenness_z y1972#c.betweenness_z y1973#c.betweenness_z y1974#c.betweenness_z y1975#c.betweenness_z ///
			y1976#c.betweenness_z y1977#c.betweenness_z y1978#c.betweenness_z y1979#c.betweenness_z y1980#c.betweenness_z ///
			y1981#c.betweenness_z y1982#c.betweenness_z y1983#c.betweenness_z y1984#c.betweenness_z y1985#c.betweenness_z ///
			y1986#c.betweenness_z y1987#c.betweenness_z alt_1988 y1989#c.betweenness_z y1990#c.betweenness_z ///
			y1991#c.betweenness_z y1992#c.betweenness_z y1993#c.betweenness_z y1994#c.betweenness_z y1995#c.betweenness_z ///
			y1996#c.betweenness_z y1997#c.betweenness_z y1998#c.betweenness_z y1999#c.betweenness_z y2000#c.betweenness_z ///
			y2001#c.betweenness_z y2002#c.betweenness_z y2003#c.betweenness_z y2004#c.betweenness_z y2005#c.betweenness_z ///
			y2006#c.betweenness_z y2007#c.betweenness_z ///
			i.year i.district_code  i.year##i.depa ,  cluster(district_code)

eststo
estadd scalar obs_shown = `obs_shown'



capture cd "$tables"
esttab est* using "tableA20.tex", replace nonumber label style(tex) ///
		cells(b(star fmt(%12.4fc)) se(par fmt(%12.4fc))) ///
		 stats( obs_shown r2 , fmt( %12.0fc %12.4fc ) ///
            label("\hline Observations (1983-1995)"  ///
			  "\$R^2\$" )) ///
		keep( 1.y1983#c.betweenness_z 1.y1984#c.betweenness_z 1.y1985#c.betweenness_z ///
			1.y1986#c.betweenness_z 1.y1987#c.betweenness_z   1.y1989#c.betweenness_z ///
			1.y1990#c.betweenness_z 1.y1991#c.betweenness_z 1.y1992#c.betweenness_z 1.y1993#c.betweenness_z ///
			1.y1994#c.betweenness_z 1.y1995#c.betweenness_z ) ///	
		order(1.y1983#c.betweenness_z 1.y1984#c.betweenness_z 1.y1985#c.betweenness_z ///
			1.y1986#c.betweenness_z 1.y1987#c.betweenness_z   1.y1989#c.betweenness_z ///
			1.y1990#c.betweenness_z 1.y1991#c.betweenness_z 1.y1992#c.betweenness_z 1.y1993#c.betweenness_z ///
			1.y1994#c.betweenness_z 1.y1995#c.betweenness_z) ///
		varlabels(1.y1983#c.betweenness_z "Year 1983 \$\times\$ Average Betweenness" 1.y1984#c.betweenness_z "Year 1984 \$\times\$ Average Betweenness" ///
		1.y1985#c.betweenness_z "Year 1985 \$\times\$ Average Betweenness" ///
			1.y1986#c.betweenness_z "Year 1986 \$\times\$ Average Betweenness" 1.y1987#c.betweenness_z  "Year 1987 \$\times\$ Average Betweenness" ///
			1.y1989#c.betweenness_z "Year 1989 \$\times\$ Average Betweenness" ///
			1.y1990#c.betweenness_z "Year 1990 \$\times\$ Average Betweenness" 1.y1991#c.betweenness_z "Year 1991 \$\times\$ Average Betweenness" ///
			1.y1992#c.betweenness_z "Year 1992 \$\times\$ Average Betweenness" 1.y1993#c.betweenness_z "Year 1993 \$\times\$ Average Betweenness" ///
			1.y1994#c.betweenness_z "Year 1994 \$\times\$ Average Betweenness" 1.y1995#c.betweenness_z "Year 1995 \$\times\$ Average Betweenness") ///	
		collabels(, none) eqlabels(, none) starlevels(* 0.1 ** 0.05 *** 0.01) ///
		prehead(\begin{tabular}{l*{@M}{c}} \hline ///
			&  \multicolumn{1}{c}{Land Occupations (log)} &  \multicolumn{1}{c}{Land Occupations (binary)}   &  \multicolumn{1}{c}{Squatter Peasants (log)}  ///
			&  \multicolumn{1}{c}{Occupied Hectares (log)}  \\ ///
			& (1)  & (2) & (3) & (4)              \\)  ///
		posthead(\hline ) mlabels(, none) ///
		prefoot( ) postfoot(\hline \end{tabular})
